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INTRODUCTION 


The  primary  purpose  of  linear  regression  is  the  prediction  of  Y 


from  a given  value  of  X by  means  of  the  regression  equation  Y=BX  or 


and  the  estimation  of  confidence  limits  for  true  Y 


Regression  analysis  for  two  variables  may  be  made  on  a weighted  or 


unweighted  basis  with  the  line  passing  through  the  origin  or  through  the 


means 


NON-CONSTANT  VARIANCE  OF  Y AND  WEIGHTING 


If  Y is  the  dependent  variable,  the  linear  regression  of  Y on 


X is  measured  without  ervo 


measurement  of  Y has  the  same  variance  over  the  whole  range  of  X 


When  this  is  not  the  case  deviations  from  the  regression  line  will  not 


be  normally  distributed  and  the  prediction  of  Y from  X will  be 
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When  the  variance  requirement  is  not  met,  a transformation  of  Y 
to  logarithms  or  square  root,  etc.,  may  be  tried  whereby  the  variances 
at  the  lower  end  of  the  scale  of  X will  be  found  to  be  non-signifi- 
cantly  different  from  those  at  the  upper.  If  repeat  measurements  of 
Y at  one  value  of  X are  available,  the  reciprocal  of  the  variance 
derived  from  such  a cluster  may  be  used  as  a weight  in  the  analysis. 

If  variable  Y comprises  both  single  measurements  and  averaged 
measurements , the  number  of  measurements  that  were  averaged  should  be 
used  as  a weight.  If  both  the  variances  of  Y and  number  of  measurements 
per  averaged  Y are  known  the  weights  would  be  the  number  of  measurements 


■'i 

L.  * J 

. :♦< 

ft 


of  Y at  X^  divided  by  the  variance  of  Y at  X^  . If  the  variance 

_2 

of  Y at  X^  is  a function  of  X then  X^  may  be  used  as  a weight. 

PROGRAM  EXTENSIONS 

Program  REGRESS  (S1281)  was  developed  and  written  by  J.  Skory  and 
Z.  Ungar  at  the  United  Aircraft  Corporate  Systems  Center  for  the  IBM 
1800.  It  has  been  adapted  to  the  UNIVAC  1108  with  the  following 
modifications  and  additions: 

(l)  the  maximum  input  sample  size  has  been  increased  to  2000 
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with  k iuternled  i'rom  1 to  u for  a sample  of  n 

(3)  Confidence  limits  for  true  Y are  calculated  at  the 

experimental  X?  and  post  experimental  X^  . 

(4)  Effective  N at  Xj  is  calculated  for  use  in  finding 

tolerance  limits  for  Y . Effective  N in  a value  such 

that  average  Y at  a given  X is  predicted  as  accurately 
from  the  point  on  the  regression  line  as  if  N observations 
had  been  actually  made  «t  that  specified  X . 

(5)  Linear  calibration  formulas  are  given  for  prediction  of  X 
from  Y with  confidence  limits  for  true  X 

(6)  Formulas  are  given  for  the  statistics  necessary  for  the 
comparison  of  2 regression  lines.  Test:  are  provided 
for  homogeneity  of  variances  about  t.he  lines,  parallelism 
and  coincidence. 
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(The  last  three  items  are  calculated  by  the  user  after  the 
• program  has  been  run.) 

INPUT  INFORMATION 

The  program  reads  all  the  following  input  data  from  the  card 
reader  (unit  3) : (For  I/O  formats  see  Appendix  A) . 

(1)  n the  number  of  data  points  in  the  input  sample. 

(2)  IH  an  indicator,  = 1 for  regression  through  the  origin 

(Y  = BX)  or  = 2 for  regression  through  the  means  (Y  =A+BX). 

(3)  IW  an  indicator,  = 1 if  all  weights  = 1 , or  = 2 if 

-2 

Wf  = Xf  , or  = 3 if  weights  are  given. 

. (4)  t Student's  t for  calculating  confidence  limits  fbr 

true  Y . 

(5)  X^...Xn  the  independent  variable  at  n points  equally 
or  unequally  spaced  along  the  X axis  and  assumed  to  be 
controlled  by  the  experimenter  without  error  (or  only  with 
negligible  error). 

(6)  Y^...Yn  the  dependent  variable.  Y^  is  an  associate  of 

Xj_  and  if  Y is  an  individual  single  measurement  Var(Y^)  = 
Var(Y2)  = ...  = Var(Yn). 

(7)  Wi...Wn  weights  for  Y^  given  as  input  only  if  IH  = 3. 

(8)  ?!• ••?!()  ten  points  of  X chosen  post  experimentally  at  which 
it  is  desired  to  predict  Y and  estimate  its  error  as  an 
exploratory  step  for  future  experimentation. 

I 
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(9)  Q^...Q10  ten  weights  associated  with  P.^  . 

COMPUTER  ANALYSIS 

Program  REGRESS  first  sorts  the  points  of  the  independent  variable 
X into  ascending  order  using  a modified  Shell  sort  embodied  in  S0RT,  a 
Fortran  subroutine  (NUSL  Program  1525)  programmed  by  Lesick  and  Shores 
and  modified  to  retain  the  and  associated  with  each  X^  . This 

is  done  in  order  to  obtain  the  within  groups  sums  of  squares  for  use  in 
a one-way  analysis  of  variance. 

The  program  then  evaluates  the  equation  Y=A+BX  or  Y=BX  where 
A and  B are  weighted  least  squares  estimates  determined  by 
minimizing  the  sums  £ w:  (Yl-fy 1 or  £ Wi(Z0-  Zif  where  the 
are  translates  of  Y^  . The  program  thus  performs  a weighted  linear 
regression  analysis  through  the  means  or  through  the  origin. 

The  deviations  of  observed  Y from  the  regression  line  are  tested 
for  normality  using  WTEST  , a Fortran  function  (NUSL  Program  SI64O) , 
(References  2 and  3). 

The  program  writes  all  output  data  on  the  line  printer  (Unit  4) s 

SUMWX  and  SUMWY  sum  of  the  weighted  X^  and  Y^ 

SUMW  sum  of  the  weights 

AVERX,  AVERY  weighted  average  of  X^  and  Y^ 


SUMDX  SUMDY  SUMXY  weighted  sum  of  squares  and  products 

* about  the  mean  or  origin 

1 B the  regression  coefficient 
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the  intercept 

predicted  Y vector 

predicted  exploratory  Y vector 

vector  of  deviations  of  observed  Y from 

the  regression  line 

correlation  coefficient  squared 

standard  error  of  estimate 

a variance  ratio  testing  if  B=R=0 

F(l,n-1),  through  the  origin,  is  evaluated 

F(l,n-2),  through  the  means,  is  evaluated 

standard  error  of  B 

standard  error  of  A 

another  variance  ratio  testing  if 

A=0  (for  case  of  through  the  means). 

F(l,n-2)  is  evaulated  with  1 and  (n-2) 
degrees  of  freedom 

standard  error  of  each  predicted  point  on 
the  line  at 

standard  error  of  each  predicted  exploratory 
(postulated)  point  on  the  line  at  post 
experimenta  L Xi (-P^) 
standard  error  of  each  weighted  newly 
observed  Y or  (Y^-Y^) 
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Sy  . 
1ex 


standard  error  of  each  weighted  newly 


VARY 

SSB 

SSDEV 

U 

CL 

N 


observed  Y. 


ei 


or 


(Yei-Yei) 


Student's  t for  each  deviation  from 
regression: 

t(n-l)  for  regression  through  the  origin 
t(n-2)  for  regression  through  the  means 
variance  of  Y within  groups  at  one  X 
sum  of  squares  due  to  B 
sum  of  squares  of  deviations  from  the  line 
normal  distribution  test  statistic 
(applicable  only  when  all  are  equal) 
confidence  limits  for  true  Y at 
experimental  and  post  experimental  X^ 
effective  N 

The  mathematical  theory  and  derivation  is  presented,  for  example,  in 
K.A.  Brownlee  (i960).  In  the  calculations,  many  of  the  statistics  are 
evaluated  by  alternate  algebraically  equivalent  methods. 


POST  EXPERIMENTAL  ESTIMATION 

After  an  experiment  is  completed  one  may  wish  to  estimate  Y at 
3ome  points  of  X(=P)  not  already  tried  in  the  experiment  and  obtain 
confidence  limits  for  these  estimates.  In  addition,  one  may  wish  to 
know  how  many  repeat  measurements,  Q , are  required  at  a specified 
value  of  X(=P)  so  that  the  error  of  a weighted  newly  observed  Y is 
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less  than  0 at  95%  confidence.  This  may  be  done  by  making 
and  varying  Q1,  Q2,  widely  and  plotting  the  error  against  Q. 

POST  COMPUTER-ANALYSIS  CALCULATIONS 
From  the  computer  output  further  calculations  may  be  made  by  desk 
calculator  to  obtain  the  following  five  types  of  estimated  information. 

(1)  Tolerance  Limits 

Using  effective  N in  the  manner  of  (4)  (page  498)  the  tolerance 
limits  for  Y may  be  estimated.  The  following  is  the  type  of  statement 
that  would  result;  There  is  95%  confidence  that,  at  the  X specified, 

99 % of  the  population  of  individual  measurements  (observations)  of  Y , 
in  future  experimentation,  will  fall  between  the  limits  calculated. 

(a)  Tolerance  limits  at  experimental  equals 

Y^ru^S  (Using  N±) 

(b)  Tolerance  limits  at  post  experimental  X^ 

Yei±(ru)eiS  (Using  Nei) 

(c)  Tolerance  limits  for  W-_  newly  observed  Y at  experimental  X^ 

Yi+(ru)1S  (Using  Nj.) 

(d)  Tolerance  limits  for  newly  observed  Y at  post  experimental  X^ 

Yei±(ru)eis  (Using  Nei) 

m 

(2)  Linear  Calibration:  Prediction  of  X from  Y 

A.  Sample  Unweighted:  Line  Through  the  Means 

(l)  Confidence  limits  for  true  X at  a weighted  newly 
observed  Y 


8 


NUSL  Tech  Memo 
2070-451-69 


X^±tF'"| 


w 


where  t = Student's  t for  (n-2)  degrees  of  freedom 
F comes  from  test  of  significance  for  B . 

If  F<1  results  may  not  be  meaningful. 

W is  the  number  of  Y's  averaged  to  give  the 

new  observation  Y . 

„ upper  limit-lower  limit 

(2)  Standard  error  of  X = 

2t 

B.  Sample  Unweighted:  Line  Through  the  Origin 

(l)  Confidence  limits  for  true  X at  a weighted  newly 
observed  Y n 


W B1'  / 

with  t at  (n-l)  degrees  of  freedom. 


. . „ upper  limit-lower  limit 

(2)  Standard  error  of  X = — 

2t 

C.  Sample  Weighted:  Line  Through  the  Means 

(1)  Confidence  limits  for  true  X at  a weighted  newly 
observed  Y are  obtained  by  solving  for  X in  the  quadratic 

(Y  - Y - B ( X • - X ) ) * 1 = ^ (( $ ^ ) t K (*C - * ) : ^ - ( X • - X )‘J 

_ upper  limit-lower  limit 

(2)  Standard  error  of  X - 

2t 
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D.  Sample  Weighted: 

Line  Through  the  Origin 

(l)  Confidence  limits  for  true  X at  a weighted  newly 

observed  Y 
X 

Y 

B+tS&i)* 
W n 

. x „ „ upper  limit-lower  limit 

(2)  Standard  error  of  X = 

(3)  Comparison  of  Two  Regression  Lines  (W^=l) 

Given  the  estimates:  , 

Line  1:  Xx  Yx  Bx 

-1  4 £(V-Y,)2 

Line  2:  X2  Y2  B2 

n2  s2  f2  I (X-xzf  I(Y-Yi)z 

A.  Test  If  S^  = S2 

• 

Make  2-sided 

2 

F test  with  larger  S in  the  numerator . 

If  S^  is  larger: 

S^ 

1 < 

F = with  (fx,  f2) 

S2 

B.  Test  If  B^  = B2 

(A  Test  of  Parallelism) 

If  S^  is  not  significantly  different  from  S2  they  are 

pooled  and 

S = 

fl  + f2 

Then 

B'Z?L - 

1 1 'i 

Tv.-V*  + 
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Test  If  the  Two  Lines  Are  Coincidental  (if  X^-X2i  Yy-Y25 
given  B-^=  E^) 

Pooled  estimate  of  B 

= b 

flXi-Wi-ftXi-Xj1 

£ 

Pooled  variance  of  estimate  = 5 = 

ft  Y-w+ tn-Y,)z - B'ftix-zr+Fotrir) 


r>,  rii 


+ nxrx,r 


(4)  Analysis  of  Variance 

The  computer  output  contains  the  necessary  information  to  perform  a 

one-way  analysis  of  variance. 

(5)  Simultaneous  Tolerance  Intervals 

These  intervals  are  calculated  as  Ye^  + c**Sye^  with  c**  tabled 

and  its  use  explained  in  4.7.4  of  Reference  5. 


JOHN  SKORY 
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APPERDBC  A 


MOMQICLAIUFE  TABLE 


Fortran  Variable  Mane  4 Type 
X real  array 
T real  array 
W real  array 

P real  array 

<4  real  array 
H integer 

IH  integer 

IW  integer 
STBDRT  real 

TITLE  real  array 

SUM*  double  precision 


Meaning 

The  independent  variable  of  the  input. 

The  dependent  variable  of  the  input. 

The  vector  of  weights  associated  with 
the  independent  variable  X . 

The  vector  of  post-experimental  points 
of  the  independent  variable  X . 

The  vector  of  weights  associated  with  P 

The  number  of  observations  in  the  input 
variables  X and  T . 

An  indicator  for  regression  through  the 
means  or  through  the  origin. 

An  indicator  for  the  weights  W . 

Student's  t for  calculation  of  confident 
levels. 

An  array  used  to  store  the  Hollerith 
characters  of  the  title  card. 

The  sun  of  the  weights  W . 


. 


SUMIX 

double  precision 

The  weighted  sum  of 

X . 

SUM02 

double  precision 

The  weighted  sum  of 

X2  . 

stntnr 

double  precision 

The  weighted  sum  of 

I . 
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RDMEBCLAIGRE  TABLE  (Cont'd) 

Fortran  Variable  Bane  l»  Type  Meaning 


S0MTC2 

double  precision 

The 

SUMNXY 

double  precision 

The 

AVERX 

double  precision 

The 

AVERY 

double  precision 

The 

SUMDX 

double  precision 

The 

of 

SUMDY 

double  precision 

The 

of 

SOMXY 

double  precision 

The 

of 

B double  precisian 

The 

A double  precision 

The 

THAT 

real  array 

The 

RECRY 

real  array 

The 

weighted,  sub  of  I . 

weighted  sum  of  XT  . 

weighted  mean  of  X 

weighted  mean  of  T . 

weighted  corrected  sum  of  squares 
X . 

weighted  corrected  sum  of  squares 
Y . 

weighted  corrected  sum  of  products 
X and  Y . 

regression  coefficient. 

regression  intercept. 

predicted  Y of  regression. 

deviations  of  Y from  the  regression 


SUMRS  double  precision 

RSQP  double  precision 

SSB  double  precision 
8 double  precision 
AJfl  double  precision 
SYHAT  real  array 

SY  real  array 


The  weighted  sum  of  squares  of  the 
deviations. 

The  multiple  and  simple  correlation 
coefficient  squared. 

The  weighted  sum  of  squares  due  to  B . 
The  standard  error  of  estimate. 

Multiplier  for  intercept  variance. 

The  standard  error  of  each  predicted 
point. 

The  standard  error  of  each  predicted 
point  assuming  these  points  are  newly 
observed  Y with  weights  *»  W . 


NOMENCLATURE  TABLE  (Coat'd) 
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fortran  Variable  Hamc  t Type 
T real  array 

CLYIH  real  array 

CLY1L  real  array 

CLY2H  real  array 

GLY2L  real  array 

XHAT  real  array 
HI  real  array 

YSUBB  real  array 


Student's  t for  each  deviation 

from  the  regression  line. 

The  upper  confidence  limits  for  true  Y 
at  each  experimental  . 

The  lover  confidence  limits  for  true  Y 
at  each  experimental  ^ . 

The  upper  confidence  limits  for  true  Y 
at  each  experimental  X^  assuming  these 
points  are  newly  observed  Y with 
weights  » W . 

The  lower  confidence  limits  for  true  Y 
at  each  experimental  X^  assuming  these 
points  are  newly  observed  Y with 
weights  * W . 

Effective  M for  tolerance  limits  at 
each  experimental  Xi  . 

Effective  M for  tolerance  limits  at 
each  experimental  X.  assuming  these 
points  are  newly  observed  Y with  weights 
* W . 

The  predicted  exploratory  Y after 
regression. 


SYHE  real  array 


The  standard  error  of  each  predicted 
exploratory  point. 


SEYE  real  array 

GLYE1H  real  array 
CLYX1L  real  array 


The  standard  error  of  each  predicted 
exploratory  point  assuming  these  points 
are  newly  observed  Ye  with  weights  ■ Q 

The  upper  confidence  limits  for  true  Y 
at  each  post-experimental  Xei  . 

The  lower  confidence  limits  for  true 
at  each  post-experimental  Xei  ‘ 
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NOMENCLATURE  TABLE  (Cant'd) 


Fortran  Variable  Name  fc 
CLYE2H  real  array 


CLYE2L  real  array 


BEHAT  real  array 


HE  real  array 


F12B  double  precision 
SB  double  precision 


SA  double  precision 


U real 


StM!  integer 


SUMSQ  double  precision 
BGB0UP  integer 
SSY  double  precision 
TEAR  double  precision 


The  upper  confidence  Units  for  true  T 
at  each  post-experimental  XQi  assuming 

these  points  are  newly  observed  Ye  with 

weights  « Q . 

The  lower  confidence  Units  for  true  T 
at  each  post-experimental  Xe^  assening 

these  points  are  newly  observed  Te  with 
weights  - Q,  . 

Effective  H for  tolerance  Units  for 
each  post-experimental  Xei  * 

Effective  H for  tolerance  Units  for 
each  post-experimental  Xei  assuming 

these  points  are  new  Y0  with  weights 

- Q . 

The  F test  for  B * B ■ 0 . 

The  standard  error  of  the  regression 
coefficient . 

The  standard  error  of  the  regression 
intercept. 

The  F test  for  A = 0 . 

W-test  for  nomaUty  of  the  deviations. 
The  degrees  of  freedom  within  groups. 

The  sun  of  squares  within  groups. 

The  size  of  each  group. 

The  variance  of  each  group. 

The  mean  of  each  group. 
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APPENDIX  B 
I/O  FORMATS 

The  program  reads  (TITLE(l),  1 = 1,14)  from  the  first  card  in  13A6, 
AZ  format. 

The  program  reads  N,IHfIW,  and  STUDNT  from  the  second  oard: 


Variable 

Columns 

Format 

N 

1-5 

15 

IH 

6-10 

15 

IW 

11-15 

15 

STUDNT 

16-30 

E15.8 

The  program  then  reads  the  input  variables  X,Y,W  (only  if  IW  s 3), 

t 

P,  and  Q separately  from  cards  in  5E15.8  format. 

The  program  writes  all  output  data  on  the  line  printer  in  7E17.8 
format. 
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000001 

C 1 

000002 

v» 

000000 

l 

000004 

i 

OOuOOb 

c 

00000b 

000007 

J. 

000008 

<L 

000009 

000010 

oouoii 

000012 

1 

000013 

c 

000014 

J 

00001b 

1UUU 

000016 

000017 

000018 

000019 

000020 

000021 

000022 

000023 

10 

000024 

000026 

Id 

00o02o 

Uhijnt.i> J”— t'r-uOA/i,-,  H(Ji\  1 Wo— VAiti AoLL  wEIGHltU  RLGKt-bSION  (THROUGH 
tlic.  Ol'iulli  Hull  ihROuGH  IHE  MEANS ) 

JIf-,U4SiJh  ulUilmi  X(2(I00)»  1 (2000)  , «<2uOo)»  YHaT  (20  JO ) , 
i VSullElcGOUi  » Kl.oA  Y 1 200o ) » SYHm((200O),  SYUoOO),  T(2000)» 
d ,’UO)*  UI1I1  SEYE(lO),  SYHL(lu) 

JIhLNSIUH  UUY XHltDuO ) » c]_Y1l(2000) , CLYtlH (2000 ) , CLYElL (2000 ) » 
x CLraUuOOO) » CLT2L(2U0G)  , CLYE2H(20oO) , CLYE2L(2000) , 

c NrmTUdOO)  » NLHA  ( (cOOO  ) , NI(2000)»  wE(2000) 

KEhu  N11AT  » >,1  * nEmAT  , wt 
iNFt-GLn  sU*Mi»  U^GFlil,  (Jt-(jFC 
jOodLL  HKUUiUIOu  SUMW,  SijM„X,  SUMWX2,  SU|.wY,  SUM*Y2»  SCMXY#  SUMwXY 
x * SUMUa,  S.JMUY » SUmRS#  RSOU,  S,  SA » Sb,  At  d»  AVtRx*  AVtRY» 

t FI^a,  rl2L»»  Awl » Sc,  S j6 , SUMY » SSY#  SUhSU,  V6«K»  .uNSGb, 

J MNbUiix,  ,MjGu2»  MHS'lC,  bit  F2,  C»  bbSMSG,  SrfGKFY 

xuou  rttAJ  l J,  70u,Umj-:999)  ( T I TLL  ( 1 ) ,1  = 1,14) 
i(Lhu  ( *»»  6)  i,»  lb,  I,  iSTUONT 
rtt«j  ( O,  o)  ( A ( I ) , . = 1,1.) 
kLmU ( 3 , u ) ( 1 (I)  ,i  = l,ij) 

IF  ( Xft  .Evi.  J ) ULAJ (0,6)  (4(1), 1=1, U> 

AE.Au  ( 3 » o ) (H  ( I ) , l = 1,x0) 

AEmU  (J,L>)(u(l),l  = l,iO) 

uO  1 0 ( X6* 1 j, 1 / ) , i« 

16  uU  Xo  1-1, i. 

a(D  - 1,0 
10  v.0  I V I XKUU 
1,0  10  1 t 


000027 


It,  JO  19  I-1,  M 
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000026 

..(i)  - A(i)-**(-2> 

000029 

iy 

co.«  fiiJUL. 

000000 

1/ 

CAlu  bO'UOlA'  r».<#  1,U'  ,TrtU&t) 

000001 

A(-  ( lit  2 ) Gu  ro  27 

0000.12 

*K  1 f c.  ( 4 * 6 / ) 

000000 

GO  fu  to 

000004 

2 f 

*KiTL  ( -♦  p 7 ) 

OOOOOb 

26 

wRllt(4»ttOO)  (TIIL£U)»1=;1.14) 

000006 

oKiTt. (4r  12  J h#  IW  piTUDUT 

000007 

oKilL  (4p2.iJ 

000006 

«KiT£  (tpbw)  ( X ( 1 ) p 1 = 1 pi.) 

000009 

*R1 ft(4p6) 

000040 

-KlTL(4PbU) (Y(I) »Izl,N) 

000041 

uHl TE (4  p 24 ) 

000042 

ttkiltUpbU)  (W(I)pI=1,N) 

000040 

*Rl7t  (4p21) 

000044 

<*K1 1 b (4  p 60 ) < h-  ( 1 ) p 1=1  p 10  i 

00004b 

rfRiTE<4P2S> 

000046 

aKI  TL  (4  p 60  ) (4'  ( X ) p 1-1  p 10  ) 

000047 

bUf.,*.  = U.CJ40 

000046 

GUN..  A = O.OLMJ 

00O049 

bUl'.*A2  - U.ULU 

OOOObO 

j Ul*.  A T — 0 • v 0 

000061 

bUi-.oYt  - OpKLj 

UOOObt 

oLf-ihAT  — 0 . 0 

000060 

AVut'A  - 0 • ouO 

00U064 

— OtOuO 

OOOObb 

bOl-  UA  = U • UL'O 

QOuObb 


jUi.uf  = O.ot-u 


00u057 

000066 

00uob9 

OOUuoO 

0OOU61 

O000o2 

000063 

OOOObH 

000065 

000066 

000067 

00006b 

000069 

000070 

000071 

000072 

000073 

000074 

000075 

00b076 

000077 

000076 

000079 

000060 

000061 

00006c 

0OU063 

Q00064 

00006b 


wa  T6ch  Nmo 
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jCi'-.ai'  — 0*bub 

Ih  ( 11.  ,tv..  2 ) TO  jl 

JO  c'J  l-l#n 

bUl-i*  A-ab.MHA'Yfc  ( i > * A ( 4 ) 

bOriw a2  - So.-m»a2  ♦ « ( I ) *X ( I ) *X ( I ) 

bUp,«  Y=Sji*6  » ♦*  ( I ) * Y (4  ) 

aW.*Y2  - bOiiWYZ  ♦ Nll>*1 11)»YI1) 

SUNiwXY  £ SO. ‘WAY  ♦ k(l)*A(l)*Y(I) 
bU.'  W=SUh*+A  ( I ) 

2*  COi.T  ItNUt. 

SUiiUX  = SCM«X2 
SUpiuY  = SUrt*Y2 
SUmAY  = SU"*XY 
SO  TO  Jc 
31  OO  26  l=l».N 

SU.-ibX  = SUpioX  ♦ SUM«/(SUMW**(n  ) * W ( I ) a (X(i)~AVEHX)*«2 

SUi-ioY  = SOwUY  ♦ SUM*/  (SUMN-t* ( I ) ) * *(I)  * (Y(I)-AVEKV)**2 

SUhXY  = SUCAY  ♦ SUMN/(SuMN*t.(in  * * ( I ) * ( X ( I )-AVtKX)  * ( Y ( I )— AVErtY ) 
AVlka  S SU»*/(SufcW+*(I) ) * AVLKX  + *(l)*X(l)/( summon (1)) 
aVuKY  = SUI-.*/(S>W>«(I))  * mVEKY  f w ( 1 ) *Y  ( 1 ) / (SUM64*  ( 1)  ) 
SUmxXsSumaA** ( 1 »*X( 1 ) 

SUMM  A2  - SOi'i*  a2  ♦ Y(1)*A(1)*XIJ) 

SUi.,aY=Sopi4 Y ♦ * U » *Y  1 4 > 

aUiwYc  = Si>.*Y2  ♦ wtI)*f(l)*Y(I) 

bUi'ihxY  = bo.l*A Y ♦ *(i)«A(I)*Y(I) 

3O1  i*=So>i(**»  i 1 ) 

26  Lbr.lll.UI. 

92  jibOMl/SU.-.OX 
A=MOEhY-U*A»LKA 

20 


I 


_ .. 


.1 


00006b 


Sti-ihb  = U.O.JI) 
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000067 

OOuOtoo 

000069 

000090 

000091 

000092 

000093 

000096 

000095 


000096 

000097 

000096 

000099 

000100 

000101 

000102 

000103 


000106 

000105 

000106 

000107 


000106 

000109 

000110 

OOOlll 

000112 

000113 

000116 


jl)  1J>  l=im 
fh*i ( j);a+o*X(1) 

<€ <,KY(1>  = V ( 1 ) - tnAT(i) 
bUi-.rtS  - SU.'.HS  + MU«(Y(I)-YHA1  (I)  )**2 
10  CON I iNOt. 

t\S^U=  ( SuMU  1 -SUM  ) /6U.4C  I 
SSo=K6C»u*5Ui  .or 

a2=  ( 1 . 0U0-K5OL. ) *SUMOY/  (N-IH ) 

i=b<»K  I ( b«: ) . 

Ahllsl.O/SUMi 

IF  ( 1H  .£«*.  1 ) Akl  s 0.0U0 

00  61  1=1#N 

sYHATd)  = 5 * SORTt  AU1  ♦ IX  < 1 )-A  VtRX  > **2/SUML>X  ) 

SYU)  = b ♦ SORT!  1 */M ( 1 ) 6 AW1  ♦ <X<  D-AVERX) **2/SUMOX  ) 

1 ( * ) =KLol<Y  1 1 ) /SY  ( 1 ) 

CLYlH(l)  = I HAT ( 1 ) ♦ STODNWSYHAT ( I ) 

CLYlL(l)  = rHAT(I)  - 5T00N1*SYHAT(1) 

CLY2n ( 1 ) = ThmT(I)  ♦ bYuDHT*SYlI) 

CLYcLll)  = f HmT ( 1 ) - 5YO0H1*SY(l) 

NHaY ( 1 ) z be  / jYHaT ( I ) **2 
N1U)  = 52  / aY  ( I ) #*2 
61  COu  T INUu 

uO  60  1=1 >10 
rSuiitii)=A«-i>*f>(i) 

SYittU)  = 6 * b <KT ( AMI  ♦ (tJ ( 1 ) “AVlNX ) **« /bUMOA  ) 

bfcrt(i)  = 5 « b vRT  ( 1,/«(I)  6 A»1  t (P(lfAVLRX)**2/bUMuX  ) 

CLYEHt(l)  = YbUufc  ( 1 ) f jTUDhT *SVHL ( 1 ) 

CLrtlL(i)  = TSOjE(i)  - bTUOt,T*bYHMl)  1 


QOullS 


CLU-2ri(i)  - YbU)t(i)  ♦ STUUia*bEYEll) 
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000116 

000117 

OOUllb 

000119 

000120 

000121 

000122 

000123 

000124 

000125 

000126 

000127 

000128 

000129 

000130 

000131 

000132 

000133 

000134 

000135 

000136 

000137 

000136 

000139 

000140 

000141 

000142 

Q00143 


CLYL2L11)  = YSUiE(I)  - STUURTwbEYE ( 1 ) 
REi-iaT  ( 1 ) = S2  / SYHL  ( 1 ) *«2 
>'*£(!)  = S2  / SEYE(I)**2 
4o  COi<T  lRUE 

F12U=(H-IH)*RSQj/(1.0-RSQD) 
SB=AabCO)/omrtT (F12b) 

SA  = b * SORT  ( Atiill  ♦ AVt-RX**2/SUMDX  ) 
IF  ( IH  ,£<*.  2 ) p 12A  = A**2/SA**2 
wKlTE(4»10c)SUM,<X,,SUMWY*SUM* 

WRITE (4»  103  J A YERX»  AVERY 

•RITE (4» 104) SUMUX , SOMDY 

WRITE(4*103)SUMXY 

wRI fE(4»lOO)B,SSB 

WR1TE(4#10/)A 

rtRl  F L (4 » 1<20 ) A * J 

WRITE (4* 10O) 

WRlT£(4*50) ( YHAI ( I ) » 1=1 t N) 

«*RlTE(4f  10V) 

wHifE(4*bO>  lYSUatQ)  ,1=1,10) 
wRiTE(4,llo) 

wRl  1 t (4  » 50 ) FR£G:<Y(I)  ,1=1,N) 

«R iTfc ( 4 # 1 1 1 ) KbQo 

WRITE ( 4,  11c.)  S,SJMR5 

wR*TE (4 , llo)  lft»Fl2U 

.»R11E(4»114)BU 

wKi Il(4» 110)Sa 

U = wIEoT (RLGRY » N,  ALSu , ) 

IF ( 1H“2 ) 144 ,14  7, 147 


I 

fc  1 
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000144 

14/  *hilCl4 

1161  ltWt-12A 

000145 

14H  *Ri  1 1. (4 

11/) 

000146 

wRiliH1* 

50  J f lsl»N) 

000147 

wRl ILI4 

116) 

000146 

4RxIt(4 

50M5Vnc(l)*t=l»10) 

000149 

WRi 16(4 

119) 

OOUlbO 

«hilt(4 

50)(5V(1)»JS1»N) 

000151 

«RHt  (4 

125) 

000152 

mR116(h 

50MSEY£(1)'|S|,10) 

000155 

*K1 I£(4 

12H1/1 

000154 

*K*lt(4 

51MTU)»lsl,M) 

000155 

mRI IE (4 

124)  u 

000156 

WRilt(4 

125) 

000157 

rtRl 1£(4 

55) (ClYiH( 1)»CLY1L(i)>I=1»N) 

000156 

»'/Rilt(4 

124) 

000159 

*RiTE<4 

55) (CLY£lH( 1 ) •CLYEltf I ) * 1=1*10) 

000160 

*RiTt(4 

12o) 

000161 

WRXTE(*» 

55)  (CLY<>H(i),CLY2L(l)ti:l>N) 

000162 

«RiTE(4 

12/) 

000165 

«RilL<4 

55) (CkYc2M(l) »CuYt2u(I) *1=1*10) 

000164 

wRHt(4 

150) 

000165 

«Ri ft ( 4 

50)  (fiHAl  (l)»13i»i<j) 

000166 

*Rilc(4 

15D 

000167 

«Ri1E(4 

50)  (N£M>«T  ( 1 ) f isi  • 10 ) 

000166 

*Rl IL(H 

154) 

000169 

MHi lfc(4 

50)  (WKD'Ul'N) 

000170 

OKI  It (4 

155) 

000171 

«Rllt(4 

50)(Nt(l>»151»K) 

000172 

»Ri‘lt(4 

154) 

1 

1 
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‘ 000175 

SUi><u  — u 

000174 

bUluO  = 0.0 

000175 

NLM6KP  - 0 

000176 

5M0KPY  i o.ouo 

000177 

ObSi^bO  — U.oUO 

000176 

-« 

II 

000179 

39 

N&RuUP  = 1 

Q0016G 

iSfAKT  = 1 

000181 

3o 

I = 1 + 1 

000162 

IF  (1  .Lt.  14)  wiO  .TO  55 

000163 

IF  (NOKOUP  .61.  1)  60  10  4u 

000184 

60  TO  4i 

000165 

35 

IF  (XllaTAKl )-X(l) ) 34,  ,54 

. 000166 

HGKOUP  = N6nCUP+l 

000167 

60  TO  5o 

• 000188 

39 

IF  (N6K0UP  .Lt.  1)  GO  TO  59 

000169 

40 

bSY  = 0.000 

000190 

IBmK  - u.OOU 

000191 

bUMY  = U.Ou-J 

000192 

bOi-,0  - O.OUU 

000193 

I5T0P=i5TAKT+NGr<0UP-l 

000194 

JO  36  J=lbl4pr, IbTOP 

000195 

bbY  = bbY  ♦ SJ.W(bUMW+*(  J)  ) * 

000196 

YliArt  = 5UPi*/(bJ:iW  + *lJ)  > * YdAR 

000197 

bUf-iY  = bUM»  + N(J)*V(J) 

000196 

bU.-iw  = bliMW  + d IJ) 

000199 

50 

COl.1  INUt 

000200 

rtllhokP  = NUhGKP+1 

■ 


QOU201 


:>unbu=bJrtb^+bbY 
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000202 

iiSl  - sif/l'<i>K£)j('-l) 

000203 

aUr.i;  — jyhx+fiOHJUH-il 

000204 

oGo<*.5»»  - ♦ 5Oi"lY**2/nw><0UP 

000205 

- 5.twhPT  'f  bOrtY 

00020b 

NKirL(4»13J>  XUSTmkT)  #tiGR0uP*5UMY,5SY 

000207 

jO  10  34 

00020B 

42 

NMlH  = <4-ln 

000209 

.*NbGB  s SCMKS/Ai>S(Hm1H) 

000210 

dl  = &NULC5UMRS)  - 5N6L(SUMbO) 

000211 

UEoFul  = NMh-SuMN, 

000212 

IF  IbEGFUl  *Nt*  0)  03  10  45 

000213 

HN5UB1  = 0*000 

000214 

uO  TO  4U 

000215 

45 

MhS*dl  - Bl/OEGFbl 

00021b 

4b 

IF  (SUmw  .40.  0)  GO  TO  43 

000217 

UGShSU  =•  bO5M50  - «>*gRPY**2/  < SUMN4NUMGRP ) 

00021b 

NUMGrtP  - NUHGRP-l 

000219 

wK  i r E ( 4 • 13o  i SU  -lN , &UMS0  • NUMGRP  * bSSMSQ 

000220 

MNBOU2  = SOFiio/alMN 

000221 

IF  lOtGFbl  .LO.  0)  GO  10  4/ 

000222 

F2  = oi*SO«i*  / ( SOMbO*Ut-GFtJ  1 ) 

000223 

GO  10  44 

000224 

4o 

rtlU  1 (.  ( 4 9 15b  i 

000225 

•4M5bd2  s 0*000 

00022b 

4/ 

F2  = 0.OU0 

000227 

44 

C — 55u+S0.>'itS 

000220 

UEGPw  2 14i4iiin 

000229 

MMSgC  — C/uLuK 

000230 

..Kilt.  (4*  15/)  jSj9  bati»  F 12u9  NKIH»  bUN.Hb, 

SNbl»U»  UtGFUl.  Bit 
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000231 

000232 

000233 

000234 

00023b 

00023b 


000237 


00023b 

000239 

000240 

000241 

000242 

000243 

000244 

00024b 

000246 

000247 

000248 

000249 

000250 


000251 


000252 

00o253 

000254 

000255 

000256 

000257 

000258 

000259 


i- N.a06i  r F2»  SOiTNt  SlMSu*  MnSOiI d.t  DECFC»  C»  MNSUC 


■jO  lo  i JOU 


>V9 

/GO  FOrhAI (l3Au* A2) 

5 FOkmhT  (3J5»E15,  _•.) 

6 FOrmAT  (OE15.6) 

67  FORMAT (oOHiKEbRESSlUN  ANALYSIS  FOR  TWO  VA«IAbLtS--(THROUGH  THE  ORI 
iGlN) ) 

/ FOHmAT  (o9H1KEGReSSIUN  ANALYSIS  FOR  TWO  VARIABLES--  (ThROOGH  THE  MEA 
inS)  ) 

600  FORMAT  ( AilO  » i3A6 » A2  ) 

12  FOrMAT (//10a*4Hh  = »15»10X,SH1H  = »15»l0x»5hIW  = > 15 > 10X r 14HSTUUEN 
ll’S  I = »Ei5.o) 

23  FORMAT (//10a, 37. ilNPuT  VlCTOk  X (INDEPENDENT  VARIABLE)) 

50  FORMAT (/El/. 8) 

a FOMA  I (//10X,35NlNPDT  VeCTOk  Y (DEPENDENT  VARIADeE)) 

24  FORMAT (//1oX,24:  ilNPuT  VLC TOR  W (WEIGHTS)) 

21  FORMAT (//10X, 76M1NPUT  VlcTOr  P (PObT  EXPlRIMENTAL  tAPLOKATORY  POIN 
ITS  oF  X VArIAOLE-ALwAYS  TEN J ) 

23  i-CrMAT (//10A#4d!iPOST  EXPERIMENTAL  WEIGHTS  w ASSOCIATED  WITH  P(l)) 

52  FOrMAI i//1Oa,50  iSTUoENT'S  1 FOR  CALCDLAIiDn  oF  CONFIDENCE  LEVELS  = 
i»El4,o) 

102  FOrMA I (//// 10a# ibHSoM  Or  WX(I)  = » tlS.b » 5* » 15HSD A OF  wY(I)  = ,E15. 
Id » Da  * lCiiSOM  OF  ( 1 ) — rLlS'O) 

103  FORma  T (/  10a,  12Mi  LAN  OF  a - ,l.I5.U,24X,  12mMEAn  OF  Y = rci5,o) 

104  FORMaI  (/10a#29Hg5  OF  DEv.  plvoM  FilAn  OF  wX  - ,El3.6,/10x»29HS5  OF  D 
IwV.  From  MEmN  OF  wY  = »ci3,ti) 

103  FOkma  T ( / 10a » 24H jDM  uF  CrOSS-PROUUC 1 S = •(.l3.o> 

106  FORMA)  ( //10a » 29  iREoReSSION  COEFFICIENT  - i = ,Ei3.d*  10X,32HSSb  = S 


« 


0UU260 

00U2ol 

1)00262 

000263 

000264 

00u265 

000266 

000267 

000268 

000269 

000270 

000271 

000272 

000273 

000274 

000276 

000276 

000277 

000278 

000279 

000280 

000261 

000282 

000283 

000264 

000286 

000286 

000267 

000266 
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1UM  OF  bv.uAlM-6  Out  TO  6 = »E15.6) 

10/  F0i<MmH//10a»16hINTc.KCE*'T  = A = ifclS.8) 

120  FORMAT  (//1ua»27hEGUATI0H  = tHAf  = A + BX  = *E16.8»6H  ♦ »E15,8»2H 

x X) 

10a  FOrMAT (//10X# 18, iFREUICTuO  Y = THAT) 

109  F0.<iiAr(//10*,33hPREUICTt0  EaPLORATOKY  Y = Y SUU  E) 

110  FOhMAI C//10X,45mDEV1ATI0nS  OF  Y FROM  REGRESSION  lINE  (Y-YHAT)) 

111  FORMAT (//10a, 47HMUj.T.  And  SIMPLE  CORR.  toEFF.  SGUAREU  = R**2  = r£l 

I6.0) 

112  FORMAT (//1oa,29mSTO.  ERROR  OF  ESTIMATE  s S = >El;j.b»10Xr46HSSDEV  = 
1 SuM  SvlUAREa  OF  UEVlATIuNS  FROM  LINE  = »ElS.8) 

113  F0RMATl//lux,6HF(l,N-»Il,lh) #1X»26HAS  A TEST  IF  R = 6 = 0,  = »E15. 
la) 

114  FORMAT <//lUX#28HSTD.  ERROR  uF  b = S SUB  b = »E15.8) 

118  FORMAT  i//iUX,28.iSTD.  ERROR  OF  A = S SUB  A = ,E15.8) 

Ho  F0RKAT(//10A,6HF(1,n-iI1,1H).1X»82HAS  A TEST  IK  A (INTERCEPT)  = 0. 
1 Cr  IMA  I HuvjRESSlON  LlNL  GOt-S  THROUGH  THL  ORIGIN  = *E16.8) 

117  FORMAT (//lUA,76HSTDt  ERrOR  OF  EACH  PREDICTED  POINT  (ON  THE  LINE),( 
itliAT)  Hi  EaPERI  .ENTmL  X(D) 

116  FORMA  I (//lOAfiO/HSTu.  ERROR  OF  EACrl  PREl/jCTcD  EXPLORATORY  PuINT  (0 
IN  THE  L1NET#(YHhT  Sub  E)  AT  POST  EXPERIMENTAL  X (I)  OK  P(I)) 

119  FORMAT (//1ua»61mSTD.  ERrqR  uF  c-ACm  OBSERVED  Y»  (Y-YHAT)  AT  tXPERIM 
itNiAL  a(I) »/10X»70HASSUMlNo  THEs>E  r(I)  AFt  nE*L?  OBSERVED  INoIVIDU 
2ALS  nlTH  WlIGhTS  = <»  ( 1 ) ) 

l2a  pOnMAI  (//luA»104hbru,  ERROR  OF  EACit  UbSLl  Vc.0  Y Sjb  L OK  OF  EACH  (Y 
1 Sul.  l.  - Tn«T  Sub  E)  HT  POST  ExPER  IrtENTRL  X(l)  UR  R 1 1 ) t / 1 Ox  » 70HASS 
tUMxNG  intSi-  T(l)  ARu  nE^LY  ubSEKV-.O  INDUIJUhLS  rfllH  WtiGMTS  = 0(1 
8)  ) 


I . 


\ I 


'■} 


v 

. i 

(j 

: 


121  FOrmA)  (//lux»4ri(  (N— » lit  ill ) »iX#4lHAi»  TEST  uF  EACH  UEWIATxOn  FROM  RE 

27 


0002*1* 


000290 

000291 

000292 

000293 

000294 

00029b 

000296 

000297 

000298 

000299 

000300 

000301 

00o302 

000303 

000304 

000305 

000306 

000307 

000306 

000309 

000310 

000311 

000312 

000313 

00U314 

000315 

000316 

000317 
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loRi_ 

51  FGi\MAl  ( /l'15.w) 

122  I- o.m-mT  i//loX,30,<NOhi''ALicATluN  STATISTIC  = U = »Llb,b) 

123  FO.vHAU//10x,49mCONFIu£uCE  LIMITS  FOR  TRUE  T AT  EXPERIMlNTAL  XU)) 

124  POxMmI  l//lUA,o<i,,CONFIOEi'.CE  LIMITS  FOR  TRUE  Y AT  POST  EXPERIMENTAL 
1XU)  Ok  PUD 

126  FORM«T(//loX,49nCONFIOENCE  LIMITS  FOR  TRUE  Y AT  EXPERIMENTAL  X(I)/ 
110A»70hASSortING  THESE  YU)  ARE  NEWLY  OBSERVEO  INDIVIDUALS  WITH  WEI 
cjHTS  = W(l)) 

127  FORMAT 1//10a,62mC0NHDENCE  LIMITS  FOR  TKuE  Y AT  POST  EXPERIMENTAL 
1XU)  OK  PU)/10x»70HASSUMlNG  THESE  Y(I)  ARE  NEWLY  OBSERVED  INOIVID 
2UAlS  with  WLJ6HJS  = 0(D) 

53  FORmA) ( 3 (El 7,6»2H  ,#E17.8i2h  ))) 

130  FORMAT  (//IOaUIuLFFcCTIWE  N//10X»4HNHAT) 

131  FORMAT (//10a, lOhNHAI  SUb  E) 

13c.  F0KmAT(//10AUHk) 

133  FOhMmT (//10a»  7H.  1 SUu  E) 

13h  FOKMAl (///10Xrb4HVAKY  = VARIANCE  OF  Y WITHIN  GROUPS  AT  ONE  X,  GROU 
IP  sicL,  Y IwTAU 

135  FORMAT  (//Tlor  • Ai  X = • »t  l3,o#6X»  Ml  =• , I4,6x»  ‘SUM  Y = *»t,13.tt.bX. 

1 ’VARIANCE  = »,tl3,6,*  <•) 

130  KGr\mA  T (✓  /T1 1 * *W  lThli'i  GROUPS  DEGREES  OF  FrKlOuM  Al«0  SUM  uF  SGUARtS' 
1 *»  * JET 1EEN  GROUPS  lEGRLES  OF  FKLEuOH  AND  SUM  uh  SUUARtS'/ 

c (oa,  •*'/2o,<»Ib.l4X»El3.U*T64»  »20xt  I5*14x,tl3.6/T6<i,  •♦•) 

13/  hOKNAl Uhllct, 'aNALTSIS  OF  v ARlANCt ' /IX H31 (**•) //aX* ' SOURCE  OF  •» 

1 ' vAhlMl  10.4  • » Ua,  'UuGREfcS  uF  FKtEUOM*  »9X,  *SU"i  uF  SuuAKLS'  » luX* 

2 Mcrt'.  SjUaKL* » 13X» 'F'/lXf 131 (•*• )//Tb» ‘SLOPE  OF  TnE  LInE'» 

3 21a»  ‘1  * »20aiEi  j,o»  l4Xrt.lJ.o»  9x»E13.l//T6#  *L»E  VIATIOhS  FROM  '» 
'InL  LlLE*rlOX*Ib.cOX»Ll3.o»)4X,E13,6//U2*tGKOUP  i-.EANS  •» 
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O0u31d 

000319 

000320 

Q00S21 

000322 

000323 


0 uli.c*  »10A,Ib.20X»fcl3,8»14X*EJ3.6»9X»tl3.8//T12» 

e 'rtilrlli.  »l£tA,lb*20X*tl3.8«l‘*X*Ll3«b/lX»lol(,_,)/Tb# 

7 • |v/TAu‘J*»2cX#X3,20A,Elo.<J»14X,L13.t>) 

a3o  rOr<>^AI  (//Tii»  • Ti.tHE  ArtE  40  kEPLICATLS  U Y AT  ANY  SIN6LE  VmLUE  • 

1 'Jt  »'l 
C.NL 
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APP*  u~7.  j 


S1281  PR0GRRM  REGRESS 


START 


REAO  C TITLEC I ) .1  = 1 .14) 
READ  N.IH.IM.STUDNT 
READ  (X(I) .1=1 .N) 
REAO  (Y( I ) .1  = 1 ,N ) 


I W=3 


'READ  C PI  I J .1  = 1 . 10 ) 
READ  (0( I ) .1=1 .10) 


IEAO  (M(I). 1=1. N) 


'60  T0  ( ) . 1M 


00  1=1. N 
H( I ) = 1 .0 


:ALL  S0RT3(X.Y.M.\ 
1.N..TRUE.)  I 


NUSL  Tech  Memo  No.  2070-451-69 


U.  S.  Navy  Underwater  Sound  Laboratory 
NP24  - 38752  - 4 - 70 


Official  Photograph 


1^1, 


I H=2 


WRITE  TITLE 
FOR  RGRSS. 
THRU  MEANS 


WRITE  TITLE 
FOR  RGRSS. 
THRU  ORIGIN 


WRITE  (TITLE! I ). 1=1 .14) 
WRITE  N.IH.IW.STUONT 
WRITE  1X111,1:1 .N) 
WRITE  C YC I ) .1  = 1 .N) 
WRITE  ( W( I ) . 1 = 1 .N ) 
WRITE  C PC  I ) .1  = 1 .10) 
WRITE  ( Q( I ).  1 = 1 .10) 


SUMW= 

SUMWX: 

SUMWX: 

SUMWX2 

SUMWY: 

SUMWY2 

SUMWXY 

AVERX 

AVERY 

SUMOX 

SUMDY 

SUMXY 


O.ODO 

:0.0D0 

:0.000 

=0.000 

:0.0D0 

=0.000 

=0.000 

0.000 

0.000 

O.ODO 

O.ODO 

0.000 
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U.  S.  Navy  Underwater  Sound  Laboratory 


Official  Photograph 


DO  1=1. N 

SUHDXrSUriOX  ♦ SUMW/l  SUMW+W(  I ) ) ■ 
MCI)  * ( X ( I )-flVERX )»»2 
SUnOY=SUMOY  ♦ SUMW/( SUMW+W( I ) ) ■ 
HID  * CY(  I J-RVERY  )«»2 
SUMXYrSUMXY  ♦ SUMW/( SUMW+W(  I ) ) ■ 
W( I )■( X( I 1-RVERX )■( Y( I )-RVERY  ) 
RVERX=SUMM/( SUMW+W( I ) ) « flVERX  ♦ 
W( I )«X( I ) / C SUMW  + Wl I ) ) 
RVERY=SUMW/( SUhW+W( I ) 1 ■ RVERY  ♦ 
W(  I )«Y( I )/( SUMW  + W( I ) ) 
SUMWX=SUMWX  ♦ H( I )«X( I ) 
SUI1WX2=SUMWX2  4 H(  I)«X(  I )«X(  I ) 
SUMWY=SUMHY  ♦ W( I )«Y( I ) 
SUMHY2=SUMWY2  ♦ H(I)iYI1)iY(I) 
SUMHXY=SUMWXY  ♦ H(l)iX( I)«Y(  I ) 
SUMW=SUMW4W( I ) 


00  1 = 1. N 

SUMWX=SUMWX  ♦ W(  I ) >XC I ) 
SUMWX2=SUMWX2  4 W(  I )«X(  I )«X(  I ) 
SUMWY=SUMWY  4 H( I )«Y( I ) 
SUMWY2=SUMWY2  4 N(I)iY(I)iY(!) 
SUMWXY=SUMWXY  4 H(I)iX(1)iY(1) 
SUMW=SUMW4W( I ) 


SUM0X=SUMWX2 

SUMDY=SUMWY2 

SUMXY=SUMWXY 
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U.  S.  Navy  Underwater  Sound  Laboratory 
NP24  - 38754  - 4 - 70 


Official  Photograph 


00  1 = 1.  N 

SYHHTt I )=S«SQRT( RW1+ 

( XC I )-flVERX)««2/SUM0X) 

SY(  I )=S«SQRT( 1 */W( I )+flWl+ 

( X(  I )-flVERX )«*2/SUM0X  ) 

T ( I )=REGRY( I )/SY(  I ) 

CLY1L( I ) = YHflT ( I ) - STU0NT«SYHRT( I) 
CLYlHt I )=YHRT ( I ) ♦ STU0NT«SYHRT( I ) 
CLY2H(  I )=YHRT ( I ) ♦ STU0NT«SY(  I) 
CLY2L( I )=YHflT ( I ) - STU0NT«SY(I) 
NHflT ( I )=S2/SYHRT ( I ) 

Nil  I )=S2/SY(  I ) 
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WRITE  SUMWX.SUMWY.SUhW 
WRITE  RVERX .RVERY 
WRITE  SUMOX.SUMDY 
WRITE  SUMXY 
WRITE  B 
WRITE  R 
WRITE  fl.B 
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I 
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Official  Photograph 
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L(  ISTflRT  )=X(  I 


NGR0UP=NGR0UP+1 


NUSL  Tech  Memo  No.  2070-451-69 


U.  S.  Navy  Underwater  Sound  Laboratory  Official  Photograph 

NP24  - 38751  - 4 - 70 


I 


NGR0UP=1 


SSY=0 .000 
YBflR=0.0D0 
SUMY=O.ODO 
SUMW=O.ODO 

IST0P=lSTflRT  +NGR0UP- 1 


00  J=ISTRRT . IST0P 
SSY=SSY  ♦ SUMW/ ( SUMW+W( J ) ) « 
W(J)  « ( Y( J)-YBRR)««2 
Y8flR=SUMW/( SUMW+Wf J ) ) « YBRR  + 
W(  J )«Y( J )/( SUMW+WI J ) ) 
SUMY=SUMY+W( J )*Y( J ) 
SUMW=SUMW+W( J ) 


SUMUSQ=SUMUSQ+SSY 
SSY=SSY/( NGR0UP-1 ) 
SUMN=SUMN+NGR0UP-1 
NUMGRP=NUMGRP+ 1 
SMGRPY=SMGRPY+SUMY 
BGSMSOzBGSMSQ  + 
SUMY«*2/NGR0UP 


WRITE  X( I START  ) . 
NGR0UP . SUMYtSSY 
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